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In this paper we present an extensive study of the thermodynamic properties of the two-dimensional 
quantum Heisenberg antiferromagnet on the square lattice; the problem is tackled by the pure- 
quantum self-consistent harmonic approximation, previously applied to quantum spin systems with 
easy-plane anisotropies, modeled to fit the peculiar features of an isotropic system. Internal en- 
ergy, specific heat, correlation functions, staggered susceptibility, and correlation length are shown 
for different values of the spin, and compared with the available high-temperature expansion and 
quantum Monte Carlo results, as well as with the available experimental data. 
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I. INTRODUCTION 

The fully isotropic Heisenberg model may well be con- 
sidered the cornerstone of modern theory of magnetic 
systems; the reason for such an important role is the 
simple structure of this model's Hamiltonian, whose high 
symmetry is responsible for most of its peculiar fea- 
tures. Recent years have seen a growing interest in the 
specific case of the two-dimensional quantum Heisen- 
berg antiferromagnet (2DQHAF) on the square lattice, 
due both to its theoretically challenging properties and 
to its being the best candidate for modelling the mag- 
netic behaviour rf-fhc parent compounds of some high-T c 
superconductorsEro. 

As for the theory, the 2DQHAF cannot exhibit long 
range order (LRO) for T > because of its being a two- 
dimensional model wi±|h a continuous symmetry (Mermin 
and Wagner theoremQ); the study of the finite temper- 
ature paramagnetic phase is hence a matter of under- 
standing what kind of disorder one is dealing with, i.e., 
what kind of correlation exists amongst magnetic mo- 
ments on different sites. At T = quantum fluctuations 
make the system change from the classical-like Neel state 
to a ground state that can be rigorously provenQ to be 
ordered for 5>1; for 5=1/2 the situation is not clear 
yet, although more or less direct evidences for an or- 
dered ground state, even in this extreme quantum case, 
can be drawn from many different studies (for a review, 
see, for instance, Ref. [^), including the present one (see 
also Ref. ^|). 

The experimental activity stems from the existence of 



several real compounds whose crystal structure is such 
that the magnetic ions form parallel planes and interact 
strongly only if belonging to the same plane. As a con- 
sequence of such structure, their magnetic behaviour is 
indeed two-dimensional down to those low temperatures 
where the weak inter-plane interaction becomes relevant, 
driving the system towards a three-dimensional ordered 
phase; an antiferromagnetic Heisenberg interaction and 
a small spin value make these compounds 2DQHAF's. 
This is indeed the case of high-T c superconductors of the 
La2CuC"4 or S^CuC^C^ family (5 — 1/2), of other mag- 
nets such as the S = 1 La2Ni04 and K2MF4, and the 
5 = 5/2 Rb2MnF4. The interplane interaction in these 
compounds is several orders of magnitude smaller than 
the intra-plane one, thus offering a large temperature re- 
gion where the two-dimensional behaviour can be safely 
studied. Furthermore, their having different values of the 
spin allows a meaningful analysis of the spin dependence 
of the thermodynamic properties, which is essential if the 
interplay between thermal and quantum fluctuations is to 
be clarified. 

In order to go beyond the very first treatments (mainly, 
mean field theory and spin wave theory), that are mostly 
not satisfactory in coping with its strong nonlinearity, the 
QHAF has been tackled by the field, .theory of the quan- 
tum nonlinear a model (QNLctM)[]i§Ei Nevertheless, the 
approximations needed to reduce the actual spin model 
to the related field theory are drastic, as they usually ap- 
proximate the essential features of magnetic systems in 
solid state physics, i.e., discreteness, strong nonlinearity 
of the Hamiltonian, and appearance of the angular mo- 
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mentum operators (that is, in a classical-like language, 
sphericity of the phase space) . This inadequacy becomes 
more evident when experimental data are available and 
a quantitative comparison with theoretical predictions is 
attempted; realistic spin models and more refined meth- 
ods must then be used, and the number of fit parameters 
minimized to let the real compounds drive the overall 
comprehension of the problem. Indeed, despite the suc- 
cess in explaining some of the early experimental data, 
the QNLerM approach does not lead to a satisfactory un- 
derstanding of the problem when higher values-ef the spin 
and higher temperatures are to be consideredi9. 

This work follows a previous paperQ, and, together 
with the latter, it is our attempt to move forward study- 
ing the finite temperature properties of the 2DQHAF 
by the effective Hamiltonian method based on the 
pure-quantum self-consistent harmonic approximation 
(PQSCHA) as developed in Ref. 11 ; we report here on the 



detailed derivation of the effective spin Hamiltonian and 
more results are shown as temperature and SPW-i-|} 
are varied, .together with experimental. jdataEjO'E 
classica£30 and quantum Monte CarloEJo (MC) simu- 
lations and high-temperature expansionE3 (HTE) results. 
No best-fit procedure is involved in the comparison with 
the experimental data, as we just need to know the spin 
model, which is unambiguously defined, in the case of 
the QHAF on the square lattice, once the values of the 
spin and of the exchange integral are given. The agree- 
ment we find with the available data, together with the 
clarity of our approach, allows us to draw a comprehen- 
sive picture of the subject, including both the analysis of 
the 2DQHAF behaviour and the discussion of previous 
approaches used for the same purpose. 

In Sec. |n|we introduce the 2DQHAF model and briefly 
describe the problems encountered by alternative ap- 
proaches. In Sec. [II the effective Hamiltonian method is 
described, and the approximations involved in this spe- 
cific implementation are discussed; the effects of quan- 
tum fluctuations on the physics of the 2DQHAF are then 
analysed in terms of the quantum renormalizations in- 
troduced by the PQSCHA. Sections [fv| and [v| contain 
our results for the thermodynamic properties (internal 
energy, specific heat, correlation functions, correlation 
length, and staggered susceptibility) compared with MC 
and HTE results, as well as with the experimental data. 



Conclusions are drawn in Sec. VI 



II. THE QUANTUM HEISENBERG 
ANTIFERROMAGNET 



The two-dimensional quantum Heisenberg Antiferro- 
magnet is described by the Hamiltonian 



H 



J 



s-i-Si+d , 



(i) 



i.d 



where J is positive and the quantum spin operators Si 
satisfy \Si\ 2 = S(S+1). The index i = runs over 

the sites of a square lattice, and d represents the displace- 
ments of the 4 nearest neighbours of each site, (±1, 0) and 
(0,±1). 

The most important feature of this model is the 0(3) 
symmetry of its Hamiltonian, implying no spontaneously 
broken symmetry for T > 0: the system does not support 
LRO at finite temperature and the standard spin-wave 
theory consequently produces unphysical results. The 
existence of local alignment directions, due to the per- 
sistence of strong short-range order up to high tempera- 
tures, makes possible the definition of a properly modified 
spin-wave theoryed whose results are remarkably good. 
However, spin-wave theory remains quite an inadequate 
tool to study the thermodynamics of the 2DQHAF and 
this indeed stimulated several authors to search for alter- 
native theories. 

It is useful to see how disorder develops in the model 
described by Eq. (|l]) when temperature is switched on: 
the ground state, hereafter assumed ordered for any spin 
value, becomes unstable because of long- wavelength exci- 
tations that gradually flip the spins as one moves far from 
the chosen origin of the lattice; their energy is small and 
they do not disturb the local order of the system. They 
have a classical character in that their contribution to the 
thermodynamics is almost the same in both the classical 
and the quantum case. What makes the latter different 
is the additional local disorder introduced by the short- 
range purely quantum fluctuations, whose renormalizing 
effect decreases as temperature increases. 

The 2DQHAF model (|l]) has been studied by several 
authors in terms of the QNLerM field theory, involving 
a three-component vector field Q with the constraint 
\Q\ = 1 and spatial integrations subjected to a short- 
distance cutoff A -1 . The model depends on two pa- 
rameters: the bare spin stiffness p° s and the bare spin- 
wave velocity Co; the coupling constant turns out to be 
go = coA/p'g. Despite their names, the parameters 
and Co are not directly related with those (J and S) 
defining the 2DQHAF: this relation is indeed the weakest 
point of the QNLerM approach. 

A first link between the two models was established by 
Haldane and AffleckD under a large-.!? condition, i.e., in 
the semiclassical limit. Their mapping gives A, ps and 
Co as a -1 , JS 2 , and 2\/2 J Sa, respectively (a is the lat- 
tice spacing); it follows that go = 2v / 2S , ~ 1 , which means 
that the semiclassical limit of the spin model (S^>1) cor- 
responds to the weak coupling regime of the field theory 
(go^l)- The validity of this approach when studying 
real compounds is more than questionable, and Haldane's 
suggestion of replacing S by \/S(S + 1) does not solve 
the problem. _ 

The way Chakravarty, Halperin, and Nelson (CHN)Q 
connected the two models has greater generality. They 
used symmetry arguments to show that the long- 
wavelength physics of the QHAF must be the same of 
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that of the QNLcrM, a result that holds regardless of 
the spin value. However, they could not define the field 
theory parameters in terms of those of the spin system; 
therefore the spin stiffness and spin wave velocity are just 
phenomenological fitting parameters to be determined 
from either experiments or simulations. 

The analysis carried out by CHN on the QNLcrM 
leads to the characterization of three different regimes, 
called quantum disordered, quantum critical (QCR) and 
renormalized classical (RCR), the most striking differ- 
ence amongst them being the temperature dependence 
of the spin correlations. If go is such as to guarantee 
LRO at T = 0, the QNLcrM is in the RCR and its long- 
wavelength (i.e., low-temperature) physics is that of the 
classical model with parameters renormalized by quan- 
tum fluctuations (/5g— >ps, cq— >c) and a short-wavelength 
cutoff of order a = c/T . As far as the correlation length 
£ is concerned, the famous low-temperature two-loop re- 
sult for the renormalized classical regime is 



= C t 



(— 

\2irps 



exp (27175,5 /T) 



(2) 



where Cg is a nonuniversal coefficient. Hasenfratz and 
NiedermayerEI have subsequently calculated the leading 
correction in T/2irps to Eq. (0), and also found the co- 
efficient = ea/8, so that 



T 



Aitps 



O 



2-nps 



(3) 



although Eq. ([?]) extends the temperature region where 
Eq. (g) can be used, both of them do not work for 
intermediate or high temperature. The fact that the 
low-temperature experimental results for S — 1/2 two- 
dimensional magnetic compounds can be fitted by those 
of a 2D QNLcrM in the renormalized classical regime, has 
lead to assert that the ground state of the 2DQHAF is 
ordered also in the S =1/2 case. 

The results by CHN are of great general value and 
have been extensively used to understand the 5=1/2 
experimental data, but their approach has several sub- 
stantial drawbacks for that purpose. First of all, for pre- 
cisely given S and J of the real compound, the fundamen- 
tal parameters ps and c remain unknown; the necessary 
best fit procedure involved in their determination intro- 
duces a substantial uncertainty in the whole of the work. 
Secondly, the restriction to low temperatures cannot be 
avoided, thus making|-t|hc HTE technique, extensively 
used by Elstner et oZ.£2l, of great help and importance 
in this framework. Finally, things get worse for higher 
values of thn_spin, as has been recently pointed out by 
Elstner et alS3; further adjustments of the fit parapeters 
are necessary to reproduce the experimental datatj and 
the dependence of £ upon S cannot be analyzed because 
it is not directly addressed by the theory. Quantum MC 
simulationsE3ll3 provide quite good results in a rather 
large temperature region, but up to now they are only 
available for S = 1/2. 



As far as the different regimes of the 2DNLcrM are con- 
cerned, the RCR is the most interesting (being the one 
linked with the real S>l/2 compounds), but the QCR 
has also attracted much interest in recent years. CHN 
found that any 2DQNLcrM with an ordered ground state 
crosses over from the RCR to the QCR at sufficiently 
high temperature. This statement cannot be trivially 
extended to the 2DQHAF, as we know that the relation 
between the two models only holds for low temperature. 
Nevertheless, the lower the spin, the lower the tempera- 
ture at which such a crossover should occur, so that for 
sufficiently small S at least some signs of an intervening 
QCR-like regime could be detectedcll. As the correlation 
length of the 2DQNLcrM in the QCR is £cxa(T) = c/T, 
such a sign could be, for instance, a temperature depen- 
dence of the measured £ which becomes less pronounced 
as T increases. No such experimental evidence exists 
for pure compounds, but some doped materials show a 
similar behaviourE3; it has been argued that, by think- 
ing of the doping as causing an effective increase of the 
quantum coupling, these experimental data could tell us 
that the doped magnet is undergoing a transition of the 
same nature of that between the renormalized classical 
and the QCR of the 2DQNLcrM. This reasoning is still 
controversial and we will come back to this point at the 
end of Sec. 0. 



III. THE EFFECTIVE HAMILTONIAN 

The effective Hamiltonian approach has been success- 
fully applied in the last decade in the study of many 
different physical problems (for an extensive review see 
Ref. |23|) ; the method is based on the path-integral formal- 
ism and allows us to express the quantum statistical av- 
erage of physical observables in the form of classical-like 
phase-space integrals. Other methods in quantum statis- 
tical mechanics formally lead to this same form (see, e.g., 
Ref. [24|), but they differ from each other because of the 
different approximations used to determine the effective 
phase-space density. 

In the effective Hamiltonian method developed for a 
flat phase space, such approximation is the pure-quantum 
self-consistent harmonic approximation, whose name, 
though long, is at least self-explanatory. The PQSCHA 
does in fact separate the classical from the pure-quantum 
contribution to the thermodynamics of the system, and 
then approximates only the latter at a self-consistent har- 
monic level. This means that the classical physics is ex- 
actly described at any temperature, and so are the purely 
quantum linear effects, as the self-consistent harmonic 
approxiamtion (SCHA) only affects the pure-quantum 
nonlinear contribution. In other words, we do not re- 
nounce the exact description of the classical behaviour in 
its full nonlinearity just because we cannot deal with the 
nonlinear quantum corrections to it; we rather approxi- 
mate only the latter. This result is specially valuable in 
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the study of magnetic systems because their behaviour is 
very often characterized by long-wavelength excitations 
(such as solitons, vortices or Goldstone modes), whose 
character is indeed essentially classical. 

The procedure leading to the effective Hamiltonian in 
the magnetic casecJ can be briefly summarized as follows. 
First of all, the spin Hamiltonian must be written in a 
bosonic form through a properly chosen spin-boson trans- 
formation. Secondly, the Weyl symbolcS of the result- 
ing bosonic Hamiltonian has to be determined and the 
quantum renormalizationsliil made explicit. Finally, the 
resulting effective Hamiltonian must be put into the form 
of a classical spin Hamiltonian (by the inverse of the clas- 
sical counterpart of the spin-boson transformation used 
in the first step). The specific case of the 2DQHAF is 
described in Appendix B, where the detailed derivation 
of the effective spin Hamiltonian, that was just sketched 
in Ref. [| is reported. 

In choosing the spin-boson transformation for the 
isotropic Hamiltonian Eq. ([!]), we exclude the Villain 
transformation which is designed for models with easy- 
plane anisotropy. Furthermore, neitheip-|the Holstein- 
Primakoflfcl (HP) nor the Dyson-MaleevN (DM) trans- 
formation apparently gives an alternative, as they both 
break the symmetry of the problem. However, the bro- 
ken symmetry of the (assumed ordered) ground state of 
the 2DQHAF is restored at finite temperature by long- 
wavelength excitations whose effect, as mentioned above, 
is almost entirely taken into account already at the clas- 
sical level. As a consequence, it is almost entirely taken 
into account also by the PQSCHA, no matter what spin- 
boson transformation is used to evaluate the quantum 
rcnormalizations, as far as the bosonic effective Hamil- 
tonian is eventually put into the form of a classical spin 
Hamiltonian. This means that by using the DM or HP 
transformation in the PQSCHA framework, we break the 
0(3) symmetry, but only as far as the pure-quantum fluc- 
tuations are concerned, so that all the essential features 
due to the symmetry of the Hamiltonian are actually kept 
and the use of such transformations is consequently justi- 
fied. As for the choice between HP and DM, we employed 
the latter in Ref. |TJ; in Appendix A we show them to be 
equivalent in this context. 

An important point is the ordering problem. Whenever 
a theory prescribes a function or functional to be asso- 
ciated with a quantum operator, it should also specify 
the ordering rule to be used for that purpose, to avoid 
an unnecessary uncertainty to enter the theory. When 
dealing with spin systems, such uncertainty often mani- 
fests itself in the ambiguous definition of th e spin leng th, 
leading to an arbitrary choice between S,y/ S(S + 1), or 
others. The PQSCHA gives an unambiguous response to 
this point, by asking for the Weyl symbol of the quantum 
operator to enter the formulae for its own renormaliza- 
tion and, as we will see later, it makes the effective spin 
length 



S = S + - 



(4) 



naturally appearQ (see Appendix A). Consistently with 
its being the spin length in this formalism, S sets the 
energy scale through the combination JS 2 . We therefore 
define the reduced temperature as 



t 



7 
JS 2 



(5) 



The final result for the effective spin Hamiltonian is 



JS 2 



J2 s i' s i+d + NG(t) , 



i,d 



6 2 h 



(6) 
(7) 



with the temperature and spin dependent parameters 



V = 



V 

^2 

1 X 



A 



2St ' 



(l-7k) 1/2 (coth/ k - , 



(8) 
(9) 

(10) 



moreover, 7k = (cosfci + cosfc2)/2, N is the number of 
sites of the lattice and k = (fcijfe) is the wave vector 
in the first Brillouin zone. T> is the pure-quantum renor- 
malization coefficient, which takes the main contribution 
from the high-frequency part (short-wavelength) of the 
spin-wave spectrum, because of the appearance of the 
Langevin function. 

As for the frequencies Wk, note that in the PQSCHA 
context they just appear in the evaluation of the pure- 
quantum renormalizations. In the general PQSCHA they 
depend upon the phase-space coordinate, but such de- 
pendence, when the system has many degrees of freedom, 
would make the evaluation of the phase-space integrals 
with the effective Hamiltonian a task as time demand- 
ing as a quantum MC simulation. Therefore, in deriving 
the above formulas, we have actually introduced a low- 
coupling approximation (LCA) to make the PQSCHA 
frequencies configuration independent (see Appendix B). 
On the other hand, since the effective Hamiltonian is af- 
fected by an approximate evaluation of the renormalized 
frequencies just at a secondary level, the approximation 
is worthy and assures the final results to depend only 
weakly on the specific LCA used. 

There are in fact several ways of defining a possible 
LCA in this context. In Ref. [ll], the fundamental phase- 
space-dependent parameters appearing in the theory (in 
terms of which the frequency are defined), say A(p,q), 
were approximated in the simplest way as A(p, q) ~ A = 
A(p Q ,q ); however, in the specific case of the 2DQHAF 
the minimum energy configuration becomes unstable as 
soon as the temperature is switched on, so that a more 
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refined LCA is in fact necessary for a proper description 
of the low-temperature regime. In Ref. |ll| we suggested 
for this case that A(p, q) could be better approximated 
with its self-consistent average defined by the effective 
Hamiltonian, i.e., 

A(p,q) dpdq A(p,q) e"^« = < A(p,q)) cS , 

(11) 

where Z is the partition function relative to 7i c ff. How- 
ever, if an analytical expression for such averages is not 
available, this approximation does not actually lighten 
the burden of the numerical work. The best one can do 
analytically is to evaluate them in the framework of the 
classical SCHA. In Appendix B we show that this leads 
to two coupled equations that give us all the ingredients 
we need to actually use the effective Hamiltonian, 

. k = 4« 2 (l- 7 ^; (12) 

n 2 = l- 1 -(V + V cl )=9 2 -^, (13) 

where T> c \ — t/2n 2 represents the contribution to the fre- 
quency renormalization due to the classical part of the 
fluctuations; at variance with the pure-quantum coeffi- 
cient 2?, which is a decreasing function of temperature, 
I'd rises with t. As for the solution of the above equa- 
tions, K 2 (t) is real and positive for t<9 4 , but becomes 
unphysically complex for t > 9 A . From the explicit solu- 
tion 

«> = i[fl» + (ff*-t)i] (w) 

one sees the instability to be originated just by the clas- 
sical contribution to k 2 . This instability is typical of the 
SCHA and it is now shared by the PQSCHA because of 
the particular type of LCA chosen in order to optimize 
the low-temperature results. A closer look to Eq. Jl3] ) 
shows that what causes the instability is the contribution 
to the frequency renormalization due to spin waves with 
long wavelength. On the other hand, linear excitations 
with A>2£ are in fact unphysical as they do not exist in 
a system with no LRO; in their stead the system devel- 
ops nonlinear excitations which are fundamental in de- 
termining the thermodynamics of the system, but whose 
contribution to the frequency renormalization is seen to 
be negligible by the same arguments showing them to be 
responsible for the vanishing of the magnetization. 

Although the instability only affects the evaluation of 
the renormalized frequencies, we want to devise a reason- 
able way to treat it, in order to optimize the description 
of the temperature region where the above mentioned 
nonlinear excitations become relevant in the thermody- 
namics of the system, which is in fact the same region 
where the instability occurs. The unphysical contribu- 
tion in Eq. ( [l3|) is then pulled up by inserting a cutoff 



|k| > 7r/£, over the AFM Brillouin zone, in the evalua- 
tion of T> c \. Such a cutoff is obviously t dependent, be- 
ing £ = and it is relevant in the temperature region 
where quantum nonlinear excitations are most important 
in the system. As we will see in Sec. [y| such region seems 
to coincide with that where other authors devised anoma- 
lous behaviours, ascribed to the occurrence of quantum 
criticality in the theoretical worksEil, or to doping effects 
in the experimental onesE3. 

Coming back to the effective Hamiltonian ([]]), we see 
the quantum effects to cause (?) the appearance of the 
spin length S; (ii) the renormalization of the exchange 
integral J by the factor 9 A (t) < 1; (Hi) the introduction 
of a uniform term Q{t). 

As far as Q(t) is concerned, although it does not 
play any role in calculating thermal averages, it contains 
the essential logarithmic term that transforms the spin- 
wave contribution to the free energy from classical to 
quantumEil. At t = 0, G(0) = and the zero-temperature 
energy per site is — 26* 4 (0), higher than the classical value 
—2, because of the zero-point quantum fluctuations. We 
have already commented on the appearance of 5 as a 
consequence of a precise ordering prescription; we just 
add here that this is a genuine quantum effect, being the 
noncommutativity of quantum operators the only reason 
for an ordering rule to be necessary. 
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FIG. 1. Renormalization parameter 8 4 vs t, for (from the 
top curve) 5" = 5/2, 1, 1/2 (solid lines), and S = 0.4, 0.3, 0.2, 
0.1 (dotted lines). 

Let us now examine the exchange- integral renormaliza- 
tion embodied in the factor 9 , the parameter 9 2 depends 
on both S and t, it is well defined as far as a physical so- 
lution for k 2 is available; it attains its minimum value 
for t = 0, where 9 2 — k 2 coincides with the one-loop 
correction to the spin-wave velocity. For increasing S or 
t, 9 2 increases, going asymptotically to 1 in the classical 
(S— >oo) or high-temperature (t— >oo) limit. The instabil- 
ity value 9 2 = is not reached for physical values of the 
spin, being 9 2 > l-^S)- 1 > for S > 0. 
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The essential information we get from Eq. (||) is that 
the 2DQHAF at an actual temperature t behaves as its 
classical counterpart at an effective temperature 

tcS = eHt) ' (15) 

or, in other terms, that the energy scale is renormal- 
ized by a temperature dependent factor 4 (t). In Figs, [j] 
and ^| we show 9 A and t c g as functions of t, for different 
spins, including some S < 1/2 unphysical values (dotted 
lines). As S decreases the difference between t and t e g 
becomes more and more pronounced, because of 9 2 get- 
ting smaller, indicating larger quantum fluctuations in 
the system. 



d N s O 



-/3H B ff _ 



= (0) eff 



(16) 




t 

FIG. 2. Effective classical temperature t e s = t/9 4 (t) vs t, 
for (from the bottom curve) S = 5/2, 1, 1/2 (solid lines) and 
S = 0.4, 0.3, 0.2, 0.1 (dotted lines). 

Our theory is quantitatively meaningful just as far as 
the renormalization coefficient T> is small enough to jus- 
tify the self-consistent harmonic treatment of the pure- 
quantum effects; although this is not obviously the case 
for the S < 1/2 low-temperature regime, the dotted lines 
at least qualitatively suggest that no critical behaviour 
occurs, no matter how small the spin value. What we 
rather see is that the sharp dependence of t c g upon t 
brings to lower temperatures those features which are in- 
deed typical of the highly disordered, high-temperature 
regime of the classical model. In Sec. V C we will come 



back to this point, in relation with the effects of mag- 
netic doping, and we just recall here that any reasoning 
about models with S < 1/2 should just be considered as 
speculative when dealing with real magnets. 



IV. THERMODYNAMIC PROPERTIES 

Once the effective Hamiltonian has been determined, 
the thermodynamic properties can be derived from the 
partition function Z = J dpdqex.p(—/3'Hcs)', more gen- 
eral statistical averages are given by 



where O = 0({si}) is obtained by the quantum operator 
O following the same procedure used to determine the 
configurational part of Ti. e s (see Appendix B); by (• • -) e fit 
we hereafter mean the classical thermal average with the 
effective Hamiltonian, which equals the classical average 
at the effective temperature t c g. 

In order to evaluate the classical 2D phase-space inte- 
grals appearing in (• • ') e g according to Eq. ([l6|), a stan- 
dard classical MC simulation is perfectly suitable; be- 
ing the existent classical MC datallSO incomplete as far 
as the correlation functions are concerned, we have per- 
formed our own simulations (following the same proce- 
dure described in Ref. [30| with slight modifications of the 
overrelaxed moves made to account for the fully isotropic 
exchange) on a 256 x 256 square lattice, for 0.54<£<1.42. 
The dimension of the lattice ensures the data to be free 
from saturation effects, according to the L>6t; criterion, 
for all temperatures but the first two (0.54 and 0.56). 
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t 

FIG. 3. Energy per spin u = (H)/{NJS 2 ) vs t, for (from 
the bottom curve) S = oo, 5/2, 1 and-J/2; the symbeJs are 
quantum MC data for S =rU2 (circlestfl and trianglesa) and 
previous classical MC datatJ (squares). 



A. Internal energy and specific heat 

According to Eq. (|l6|), the internal energy per spin 
u= (H)/(NJS 2 ) is 



u(t) = 4 {t) u cl (t eS ) 



(17) 



where u c \ is the classical energy; from the results shown 
in Fig. |^ we see the quantum renormalizations to in- 
crease the energy and flatten the curve u(t) at all tem- 
peratures and more markedly for smaller spin. Quantum 
fluctuations are responsible for both effects, as they in- 
troduce additional and almost temperature-independent 
disorder, thus making the system more unstable, i.e., 
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with many different configurations thermodynamically 
relevant, already at low temperatures. 

Consistently with this picture, Fig. ^ shows the peak 
of the specific heat moving towards lower temperatures 
and decreasing in height, as S decreases. Our data for 
the specific heat at finite spin are obtained by numeri- 
cal derivation of u(t), as given by Eq. (Jl7|); because of 
the temperature dependence of the cutoff involved in the 
LCA (as described in Sec. Ill), they are affected by con- 
sistent numerical uncertainty and should not be consid- 
ered but qualitatively for S = 1/2 where such dependence 
is more pronounced. 
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FIG. 4. Specific heat per spin c = du/dt vs t, for (from 
the top curve) S = oo, 5/2, 1 and 1/2; circles as in Fig. ^| 
squares are our MC results. For the (dotted) 5 = 1/2 curve 
see comments in the text. 



B. Correlation functions 

For the correlation functions G(r) = (Si-Si+ r ), with 
r = (n, r2) any vector on the square lattice, we find (see 
Appendix B) 



G(r) -5 2 ^( S; . Si+r ) cff 



(18) 



where 0* = \—\T> v and the renormalization parameter 
V r reads 



J-W- 

SN V V 1 



l + 7k 



7k 



1/2 



coth/ k --!-)(l-cosk-r) ; 

J k 

(19) 



note that for nearest neighbors, r = d, we have 2?d = T>. 

Equation (^) shows that the quantum correlation 
functions can be obtained by multiplying the classical 
ones at the effective temperature t Q g by the temperature 
and spin-dependent renormalization factor 9^S 2 . The 
overall effect is that of a strong reduction of the spin 
correlations, as it is clear from Fig. |, where G*(r) = 



|G(|r|)|/5 2 is shown for r = (r, 0), at a fixed tempera- 
ture and for different values of the spin. 
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FIG. 5. Absolute value of the spin correlation function 
G*(r) = \G(r, 0)\/S 2 vs r, at t = 0.7 and for (from the top 
curve) S = oo, 5/2, 1, and 1/2. The dotted lines are just 
guides for the eye. 
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FIG. 6. Absolute value of the spin correlation function 
G*(r) = |G(r,0)|/S 2 vs r, at S = 1/2 and for (from the top 
curve) t pi0.45, 0.50, 0.60, 0.75; also reported are quantum 
MC datau£l (circles). The dotted lines are just guides for the 
eye. 

In order to exploit some available quantum MC data 
for S — 1/2, we have performed four classical MC simula- 
tions with the effective Hamiltonian, at the temperatures 
corresponding, via Eq. (|l5|), to those at which quantum 
MC |Simulations have been performed by Makivic and 
DingES. As seen in Fig. [| the agreement we find is 
very good, and gives us confidence to proceed towards 
the evaluation of the susceptibility and of the correlation 
length, both indeed deriving from the correlation func- 
tions. 
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V. FACING THE EXPERIMENTAL DATA 

In this last section we compare our results with the 
available experimental data and show that, at variance 
with what seems to emerge from the analysis based on 
the QNLctM approach, there is no substantial difference 
between the agreement we find for systems with S = 1/2, 
and those with S > 1. 

We consider two fundamental physical observables, the 
staggered susceptibility \ an( i the correlation length £, 
and we underline once more that our results do not con- 
tain free parameters and need nothing but the values of 
J and S to be compared with the experimental data. In 
the case of the staggered susceptibility, the latter come 
in arbitrary units so that a multiplicative factor must 
be determined by the standard least-squares procedure. 
We also remind that our results for £=1/2 and low tem- 
peratures (i<0.35) are obtained in a region where the 
PQSCHA, given the relatively large value of the renpr- 
malization coefficient V, touches its limit of validityEil. 
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FIG. 7. Staggered susceptibility \* = x/S 2 vs *i f° r (from 
the rightmost curve) S = oo, 5/2, 1, and, 1/2; the crosses 
(our simulations, see text) and the squaresLJ are classical MC 
data, the triangles are quantum MC datalij for S = 1/2. 

A. The staggered susceptibility 



spin- and temperature-dependent renormalization factor 
8% can be easily evaluated for any value of r. 

In Fig. [?] we show x* = xl S 2 as a function of tem- 
perature for different spin values. For S — oo we have 
reported our MC data, as obtained from both the gen- 
eral definition (|2|) (curve) and the classical expression 
X* = (I J2i Sj\ 2 )/3N (crosses); also reported are classical 
MC data from Ref. [If], as welLas quantum MC data re- 
cently obtained by Kim et alL3 for S = 1/2. Figure || 
shows.our results for S — 1/2 together with experimental 
datalid for the real compound Sr 2 Cu02Cl2; the quantum 
MC results by Makivic and Ding are also shown. The 
case S T-rHi 3 considered in Fig. g, where the experimen- 
tal datatjO for the compounds La 2 M04 and K2N1F4 are 
reported. The agreement between our theoretical curves 
and the experimental data is substantially equivalent in 
both the S — 1/2 and 5 = 1 case, strongly suggest- 
ing that the difficulties encountered by a QNLcrM-based 
analysis to draw a unique picture for different spin val- 
ues, derive from the inadequacy of the theory, rather than 
from an actual difference in the thermodynamics of the 
different compounds. 
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FIG. 8. Staggspad susceptibility, x* = x/S 2 vs t, for 
S = 1/2. Circleaij and trianglestJ are quantum MC data 
(the dotted lines are just-guides for the eye); the squares are 
neutron scattering datat3 for S^CuC^Cb- 



The staggered susceptibility for the 2DQHAF is de- 
fined as 



X=^EH ri+ " G(r) 



(20) 



and by the PQSCHA result Eq. (|l8|) we find 

s(s + 1) + a* 

L r=£0 



X 



S(S+l) + S 2 Y,(-Y 1+r2 St (*r*i + ,,,n 



(21) 



By using our classical MC data for (si-Si +r ) c g we then 
obtain the quantum results for all values of S, as the 
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FIG. 9. Staggered susceptibility x* = xl S 2 vs t, for S = 1. 
Circles and squares are neutron scattering data for La2Ni04 
(Ref. 0) and K 2 NiF 4 (Ref. [[|), respectively. The classical 
result (dash-dotted line) is also reported. 
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FIG. 10. Curves as explained in text, for (frpn the top 
curve) S-.— oo, S = 5/2, 1, and 1/2, the trianglesLj and the 
squarestfl are-quantum MC data for S = 1/2, the circles-are 
HTE resultsLJ, and the diamonds are classical MC dataLj. 



HTE results (i.e., ((S, T)r*£ c i [T/JS(S+1)] for S > 1 
and T>JS) is nothing but the high-T and high-S* limit 
of Eq. (H|). Furthermore, Eq. ( p2] ) represents the cor- 
rect expression of the link between the classical and the 
quantum discrete magnetic model, to be compared with 
the one devised by CHN for the QNLcrM in the renor- 
malized classical regime, leading from the classical result 
by Brezin and Zinn-JustintEl to Eq. (||). 

Equation ( |22| ) allows us to obtain the quantum corre- 
lation length from the classical one by a simple temper- 
ature scaling. In other terms, we see that each value £ 
corresponds to a temperature S) which is different 
for different spin values. This point of view is taken in 
Fig. 10, where we report our results for i(£, S), together 
with MC and HTE data. In the 5 = 1/2 case, we find 
that £(£,l/2)/i(£,oo) = 9 4 > 0.5 for £>10; as already 
pointed out at the end of Sec. Ill, such a small value 
of 8 4 makes the PQSCHA no more quantitatively reli- 
able. Nevertheless, the figure shows that the accuracy 
of our results for t(£,S) at S = 1/2 and £ ~ 100, as 
from the comparison with the quantum MC data, is still 
better than 20%, which makes them qualitatively mean- 
ingful even in the extreme quantum regionEJ. However, 
the inaccuracy in the "renormalized" temperature gives 
larger inaccuracy in the inverse function, £(i), due to its 
exponential behaviour. 
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B. Correlation length 

The correlation length is defined by the asymptotic 
behaviour for r = |r| — > oo of the correlation function, 
|G(r)| - e- r /« . From the PQSCHA result, Eq. @, 
since for large r 9 4 — > const, we get (si-Si 
which means, comparing the definition 
Eq. @, 

£(t) =&i(*cfr) , 



and using 



(22) 



i.e., the correlation length £ of the quantum model at a 
certain temperature t equals the classical one £ c i at the 
effective temperature t c s > t. It is easy to see-that the re- 
lation empirically extracted by Elstner et alS3 from their 



FIG. 11. Correlation length £ vs t, for S = 1/2 (left- 
most) and S = 1. The symbols are experimental data; 
for 5 = 1/2: 63 Cu NQR dataM (cpcles) and neutron scat- 
tering data ifer La2Cu04 (squares^) and for Sr2Cu02Ci2 
(up-trianglealHI) ; for S p-i 1: neutron scattering dataripr 
La2Ni04 (down-trianglesH) and for K 2 NiF 4 (diamondsB). 
The classical result (dash-dotted line) is also reported; the 
dashed lines are the low-t and high-t results of the PQSCHA 
(see text). 

The experimental data for 5 = 1/2 and S = 1 are com- 
pared with our theory in Fig. [HJ at variance with what 
happens by using the QNLcrM approach, the agreement 
does not get worse in the S = 1 case. For each spin we 
have also reported the curves (dashed lines) obtained by 







using the standard LCA (high-i curves) as well as the 
more accurate one described in Sec. [II (low-t curves) 
without cutoff and hence truncated at t = 9 4 . They are 
smoothly connected by the continuous curves, i.e., those 
obtained by using the cutoff. For each spin value, the 
temperature interval where the cutoff is relevant identi- 
fies the region where nonlinear excitations are present in 
the system, and hence suggests where a possible crossover 
towards the QCR could occur. 

In order to better understand why the QNLcrM theory 
fails in describing the 2DQHAF when used to interpret 
the available reference data for 5 > 1, the function y{t) = 
t ln£ is shown in Fig. 12 as a function of t. We remind the 
reader that the QNLcrM approach gives y q (t) = mt+n in 
the quantum case (with m and n constant in temperature 
but 5 dependent), being such expression based on the 
classical two-loop result y c \{t) = £ln(/it)+z/ (with /i and 
v constants); first of all, one must hence check whether 
or not this classical result can be safely extended to the 
magnetic system. 
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FIG. 12. The function y(t) =t\n£vst, for (from-the right- 
most curve} S = oo, 5/2, 1 and 1/2; the trianglestfl and the 
diamondaifj are quantum MC data for S = 1/2; also reported 
are neutron scattering data for La2Ni04 (circlestij) and for 
K2NIF4 (squarealj). The inset reports, together with our re- 
sult for 5 = 1/2 (line), 63 Cu NQR data^circlesH) and neutron 
scatterina-data for La2Cu04 (squarealj) and for Sr2Cu02Cl2 
(trianglestj). 

Let us concentrate on the 5 = 00 (classical) MC 
data reported in Fig. |lj (rightmost curve): The com- 
plete curve cannot be fitted with a function of the form 
y c \ and such fit is solely possible if one restricts himself 
to a limited temperature range, as in the intermediate- 
temperature region there is a clear change in the slope. 
On the other hand, this should not surprise, as iln£ = 
t ln(/zi) + v is a two-loop result and it is hence bound to 
be correct only at lowest temperatures. We now expect 
this problem to propagate to the quantum case, which is 
based on the classical result. 

Let us then look at the S = 5/2 and S = 1 curves: ac- 



cording to the QNLcrM these should be straight lines, but 
they are not; actually one can easily see a change in the 
slope at intermediate temperature, followed by a curva- 
ture inversion at lower t. As in the classical case, hence, 
good fits with the function y q can be obtained either in 
the low- or (misleadingly) in the high-temperature re- 
gion, but not on the whole temperature range. To fit the 
experimental data for 5 = 1 with a straight line is se cm-, 
ingly impossible, as pointed out by several authorslllJli3 - ll3 
in the past few years. 

Finally, if we look at the 5=1/2 case it becomes clear 
why the QNLcrM approach gave such a good agreement 
when first used to fit the experimental data. The change 
in both the slope and the curvature of t In £ is less pro- 
nounced and possibly occurs at lower temperatures, the 
lower the spin: in the 5 = 1/2 case, we find difficult to say 
whether these features are still present or not, but, if yes, 
they occur in a temperature region where the extremely 
high value of £ (« 10 4 ) makes both the experimental and 
the simulation data more difficult to be obtained. The 
experimental data, as well as our results, do actually sug- 
gest a change in the slope; on the other hand quantum 
MC data by Kim et alSM, do not give evidence of such 
change. r— ■ 

As for the very-low temperature data by Suh et aZ.E3, 
it is to be noticed that the authors explained the change 
in the slope of their curve in terms of a crossover from 
an isotropic towards an easy-plane (Kosterlitz-Thouless- 
like) behaviour. We cannot question their interpretation, 
but we underline that their conclusion is based on the 
assumption that the tln£ curve for the 2DQHAF is a 
straight line, which is, as we are seeing, at least ques- 
tionable. 



C. Magnetic doping and quantum criticality 

We have already pointed out that at intermediate tem- 
peratures there is an interval, whose width is larger the 
smaller the spin, where quantum nonlinear effects (due to 
higher order terms in the coupling) are significant; such 
interval is identified, in our theory, by the temperature 
region where the cutoff is relevant, a region which is very 
easily recognizable in Fig. O (0.45 < t < 0.9 for 5=1/2 
and 0.65 < t < 1 for 5 = 1 jand seems indeed to coincide 
with that where Chubukov and SachdevEfJ suggested the 
occurrence of the QCR in the QNLcrM. 

Experimental data for magnetically doped materials 
have been interpreted by several authors in terms of 
a possible crossover from the classical renormalized to- 
wards the QCR. Let us concentrate, in particular, on 
the data by Carretta et al!EB, obtained by a scaling 
analysis of their 63 Cu nuclear quadrupole relaxation 
(NQR) data, for the correlation length of La2Cu04 doped 
with nonmagnetic impurities, i.e., for the compound 
La2Cui_ 2; Zn a ;04. In Fig. [13| we report their data for 
x = 0.018: it is evident that the doping causes a strong 
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reduction of the spin correlation, as well as a clear flat- 
tening of ^ as a function of temperature. Such flattening 
does occur in the same intermediate temperature region 
where the QCR has been suggested to occur. 

The authors have interpreted their data in terms of 
an effective reduction of the spin stiffness, consistently 
with our theory which shows similar effects to be caused 
by an effective reduction of the spin value. The curve 
in Fig. [13] is obtained by the PQSCHA with S = 0.35, 
a value which has been empirically determined in order 
to optimize the agreement with the experimental data. 
Although the use of our approach for S < 1/2 is not fully 
justified, these results do at least qualitatively suggest 
that for S < 1/2 the magnetic model moves towards 
a regime where nonlinear effects become relevant in a 
wider temperature region (the one where £ shows a clear 
plateau) and such behaviour could be seen as a signature 
of the crossover towards a QCR. 




t 

FIG. 13. Experimental data for the correlation length, for 
the magnetically doped compound I^Cui-^Zn^Gi, from 
Ref. The leftmost lines refer to S — 0.35, with the dashed 
lines as in Fig. [ll]; for comparison the S = 1/2 curve is also 
reported (rightmost curve). 



VI. CONCLUSIONS 

In this paper we have applied the-pure- quantum self- 
consistent harmonic approximation^ (PQSCHA) to the 
study of the thermodynamics of the two-dimensional 
Heisenberg antiferromagnet on the square lattice. The 
PQSCHA allowed us to reduce the evaluation of quantum 
averages to the calculation of classical-like phase-space 
integrals. Therefore, using classical MC simulations, we 
have been able to obtain results for several quantum ther- 
modynamic quantities. What is remarkable is that these 
results are fully determined by the system's parameters, 
i.e., the spin value S and the ratio between the temper- 
ature and the exchange energy constant Tj J. 

The main effect is seen to be the temperature-depen- 



dent weakening of the effective classical exchange con- 
stant, so that an effective classical temperature natu- 
rally arises, Eq. (|l5|); this leads to a very simple rela- 
tion, Eq. (B2J), directly giving the quantum correlation 
length in terms of its classical counterpart. For other 
thermodynamic quantities, such as correlation functions 
and the staggered susceptibility, the expressions involve 
further pure-quantum renormalization factors that can 
be straightforwardly computed. Even pushing down the 
theory to the extreme quantum case, S = 1/2, we find 
agreement both with quantum MC and experimental 
data, until the renormalization parameters of the the- 
ory become too large for T/J < 0.35 making the results 
of only qualitative value. 

For higher spin, S > 1, the theory is reliable at any 
temperature, and it agrees indeed with the available ex- 
perimental data, while quantum MC simulations have 
not yet been feasible. It would be interesting to compare 
also with the announceoEj experimental results for the 
5 = 5/2 square-lattice antiferromagnet Rb2MnF4, which 
are not yet available. 
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APPENDIX A: SPIN-BOSON 
TRANSFORMATION AND WEYL ORDERING 

In this appendix we compare rihe Dyson-Maleev§l 
(DM) and the Holstein-PrimakoflEZI (HP) transforma- 
tions as far as the ordering problem is concerned. The 
DM transformation for the spin operators S f± = S x ±iS v 
and S z , in terms of bosonic operators (a\ a) is 



S+ = (2S) 1/2 a , 

§- = (2S) - * a f (2S*-a t a) 

S z = S - a f a . 



(Al) 



The transformation is canonical, as from [a, fv] — 1 the 
spin commutation relations follow, with \S\ 2 = S(S + 1). 
Their Weyl symbols are found to beO 

S+ = {2S) 1/2 a , 

S~ = {2S)~* {2S-a*a) a* , 
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S z = S-a*a , 



(A2) 



For the use of the DM transformation in the 2DQHAF 
context we refer to Ref. ||. 

The HP transformation reads instead 



5+ = {2S-a)a) 1 / 2 a 

s- = {s+y , 

S* = S-a)a , 



(A3) 



and can be rewritten in terms of phase-space operators 
q = (a t +a)/(25) 1 / 2 andp = i(tf-a) / '(25) 1 / 2 , where 5 = 
S+l/2 and with commutator [q,p] = i/S, as 









5+ 


= s 


1 - 


§- 






s z 


= s 


0- 



q 2 +f + 



1 



1/2 



(q + ip) 



q 2 + p 2 



(A4) 



We have now to determine the Weyl symbols of these 
operators; indicating with Aw the Weyl symbol of the 
operator A, the following product relatione^ holds: 



(AB)w = A w exp 



25 



(d p dq - d q d f 



B 



w 



(A5) 



Let us start with the operator S + , which has indeed the 
form of the product of two Weyl-ordered operators, 5 + = 
SAB, with 



A w = 



l-](q 2 +p 2 + 1 
4 V 5 



1/2 



= 5> (<,"+*»)» 

(A6) 



and Bw = q + ip , where the coefficients c n ari se f rom 
the expansion of the square root. By using Eq. ( |A5| ) we 
then find 



(AB) W = c " + p 2 ) n e-^* d ^« d ^l 2S {q+ip) 

n 

= J2c n (q 2 +p 2 T 



= E c « 



(q 2 + p 2 ) n ~ ^(q 2 + P 2 )"- 1 
S 



(q+ip) 



1-- A {q 2 +P 2 ) 



-,1/2 



(q+ip) , (A7) 
where terms up to the first order in 1/5 have been kept; 



hence, being S~ — (S + )* and the operator S z in Eq. ( A4) 
already Weyl ordered, the Weyl symbols of the HP oper- 
ators are 



S ± = 5 



l-^q'+p 2 ) 



1/2 



(q ± ip) , 



S z = 5 1 



q 2 + p 2 



(A8) 



^ appears as the natural 



Note that the value 5 = 5 

spin length of the theory, being" |5| 2 = (5 Z ) 2 + 5+5~ = 
S 2 , both in the DM case (|a|) and in the HP case (|a|). 

In the above derivation for_the HP case we have ne- 
glected terms of the order 1/5 2 and higher: this approx- 
imation is necessary to obtain reasonably simple expres- 
sions of S w . Although no such approximation seems to 
appear when using the DM transformation, this is in fact 
a wrong conclusion. Indeed, it is well known that, in or- 
der to apply the DM transformation, and properly take 
into account the kinematic interaction, one should use, 
rather than the simple transformed Hamiltonian 7Ydm, 
the operator PHdmP where P is the projector on the 
Hilbert subspace of the spin system. Such subspace is 
generated by the eigenstates of the operator a) a with 
positive eigenvalues n<25, i.e., by the eigenstates of the 
operator z 2 = (q 2 + p 2 )/2 with 25 equispaced positive 
eigenvalues < z 2 < (2 — 1/25). On the other hand, to 
determine the explicit form of the Weyl symbol for the 
operator PTLy^mP is an impossible task, unless one only 
keeps terms up to the first order in 1/5; this means to ap- 
proximate P with the identity operator, which is in fact 
what we have done in Sec. [jlj. Both transformations can 
hence be used in the framework of the PQSCHA and do 
actually involve the same semiclassical approximation. 



APPENDIX B: THE EFFECTIVE SPIN 
HAMILTONIAN 

The general expression given by the PQSCHA for the 
LCA effective Hamiltonian istil 

H cff = e A W(p, g )-£a k c, 2 + i£ln^A, (Bl) 
k p k /k 

where Tt(p, q) is the Weyl symbol of the original Hamil- 
tonian TL and 

ij 

^r^E^* cosk.(i-j), 

k 

D< f = JtY, a ^ cosk.(i-j), 



Ok 



N 

h 
2uj v 



£ k = coth /k ~r~ i 

Jk 
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/k = 



4 = i£ e *' (i - j) 4> 



where i and i are sites on the lattice, while A 1 - and £> ;i 
. . u y 

are the LCA approximations of the fundamental renor- 

malization parameters of the theory: 

4(p,g) = 9 Pi 9 Pj e A H(p,g)~4 , 
Bf 3 (p,q) = d qi d q3 e A H(p,q) ~ B% ; 

their dependence upon the phase-space coordinate (p, q) 
is eliminated by the LCA. 

As suggested in Ref. [ll], and for the reasons given in 
Sec. Ill, in this work we define the specific LCA to be 
used, by setting A^p^q) rs {A?-(p, <?))fg HA [and sim- 
ilarly for B?-(p,q)}, where (• • -)^§ HA is the SCHA ap- 
proximation of the classica l-lik e average (• • -) c ff with the 
effective Hamiltonian, Eq. ([B8|) below. We then have (see 
Ref. Illh 



SCHA 
off 

SCHA 



Al = (d pl d P} e A H(p,q)) 
Bl = {d q .d q .e A U{p,q))l K s 

Let us consider now the case of the spin Hamilto- 
nian (Q). We consider the lattice as subdivided into the 
usual AFM positive and negative sublattices and here- 
after use the notation (— )' = ±1 for the site i belonging 
to the former or the latter, respectively. The Weyl sym- 
bol for Ti. (using the DM transformation) is given bytil 



H 

JS 2 



4e 



i,d 

1-^ 



(1 - zf)(l - zf +d ) 

2 , „2 \ 
i+d 



(<7i<7i+d - PiPi+d) 



. z? - z 2 

-i (-)'- — (<liPi+d + Piq-i+d) 



(B2) 



where (pi, qi) are the Weyl symbols of canonical operators 
such that [qi,Pj\ = i£ -1 <Jy, and z? = (q? +p?)/2. 
It is simpler to use JS 2 as the energy unit, i.e., to apply 
the above jramework to the dimensionless Hamiltonian 
Ti(p,q)/JS 2 , so that all the relevant quantities are di- 
mensionless; in particular, j3 — *■ 1/t and h — > 1/S. 

The first step to close the self-consistent scheme de- 
scribed above is then the evaluation of e A Ti(p,q). By 
performing the transformation k — > (7r, it) — k one can 
establish the identity 



i,i+r 



r 2 n (pp) 
i,i+r 



(B3) 



where r = (7+,^); in particular, only the following 
two renormalization parameters appear in the effective 
Hamiltonian: 



D = D--'^ = D.F^ 



D' = D 



(qq) 



i.i+d 



-Dr. 



(pp) 



-d J 



(B4) 



where d = (±1,0) or (0, ±1) is a nearest-neighbour dis- 
placement. We then find, for nearest neighbours i and 
j = i + d, 



A 2 2 

e zt = zt 



e A qiqj 
e A piPj 



. D , 
<Mj + D' 
PiPj 



D' 



furthermore, 



A 2 2 

e z.z. 



z?z?+D(z?+z?) + D'(q i q i -p iPi )+D 2 +D n 



and 



e A (zf+zf)(q i q i -p i pj) = (zf+zf+^iqiq^piPj) 

+AD'(zf+z?)+8DD' , 
e A (z?-z?)(q i q i +p iPi ) = (zf-zf^q-^+pipj) . 

It is important that the imaginary part of H(p, q) does 
not contribute any renormalization term: this in fact as- 
sures the final effective spin Hamiltonian, i.e., the one 
obtained after having performed the classical version of 
the inverse of the DM transformation on Ti. e s(p,q), to 
be real. By defining 6 2 = l-V/2 with V = (D-D'), we 
find 



e A H(p,q) 1 



JS 2 



2 ^ 

i,d 



'(z i + z i+d ) + z i z i+d 



4 + zf+d 



i+d 



(<Mi+d - P\P\+d) 



(?iPi+d +Pi9i+d) 



this is the only term of 7i c ff (p, q) we need to determine a k 
and b 2 , (as configuration-independent terms do not enter 
the evaluation of (• • -);?g HA ), which are easily found to 
be 



4k 2 (1 
4 K 2 (1 



7k) 
7k) 



(B5) 



where 7k = (cosfci + cos/c2)/2 and we introduced the 
SCHA renormalization parameter 



1 



1 



2NS^ 



(i- 7k y /2 coth/ k , 



(B6) 



that can also be written as in Eq. (13); the renormalized 
frequencies Wk are hence given by Eq. (12); the dimen- 
sionless parameter /k can be written as /k = Wk/(25f). 
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With the above determinations one can express T> as 
in Eq. (^), and the procedure is completed by the self- 
consistent solution of Eqs. ([p]) and As for the first 
uniform term appearing in Eq. ( B 1 ) , one easily finds 



(B7) 



To recast e A Tt(p, q) in the form of a spin Hamiltonian, 
we scale (p, q) — > (8p, 6q) so that 



e A H(p,q) 6* 



JS 2 



tE 



(1 - z?)(l - zf +d ) 



1- -!- 



z 2 - z 2 
\i Z i Z i+d 



(qiqi+d - PiPi+d) 



(qiPi+d + Piqi+d) 



this equation has the same func tion al form of the Weyl 
symbol for the Hamiltonian, Eq. (B2), so that performing 
the inverse of the classical DM transformation we eventu- 
ally find the form reported in Eq. (^), where the change 
in the phase-space measure by the factor 9 2N due to the 
above scaling is absorbed in the effective Hamiltonian as 
an additive logarithmic term, — Nt\n6 2 . 

Finally, the PQSCHA expresses the general thermal 
average of an observable 0(p, q) as the classical-like av- 
erage with the effective Hamiltonian 

(0) = i f^d( P , q) e-^(P^^ { d) eS , 



Z J (2ir) N 

of the phase-space function 

6(p,q) = e A 0(p,q) 



(B8) 
(B9) 



where 0(p, q) is the Weyl symbol of the operator O. In 
terms of the classical spin variables Eq. (B8) becomes 
just Eq. (|l6|). In the case of the correlation functions 
G(r) = {S\-S- 1+r ) the procedure leading to the Weyl sym- 
bol SVi+i^p, q) of the operator SvSi+r and hence to the 
expression of e A Si4 +r (p, q) is obviously analogous to the 
one described for the Hamiltonian, and the final result is 
easily found to be as reported in Eq. dlgl). 
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